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ABSTRACT 

Stellar haloes and thick discs are tracer populations and make only a modest con- 
tribution to the overall gravity field. Here, we exploit the virial theorem to provide 
formulae for the ratio of the globally-averaged equatorial to vertical velocity dispersion 
of a tracer population in spherical and flattened dark matter haloes. This gives virial 
sequences of possible physical models in the plane of global anisotropy and flattening. 
The tracer may have any density distribution, although there are particularly sim- 
ple results for power-laws and exponentials. We prove the Flattening Theorem: for a 
■ spheroidally stratified tracer density with axis ratio q in equilibrium in a dark density 

O [ potential with axis ratio g, the ratio of globally averaged equatorial to vertical velocity 

dispersion depends only on the combination q/g. As the stellar halo density and ve- 
locity dispersion of the Milky Way are accessible to observations, this provides a new 
method for measuring the flattening of the dark matter distribution in our Galaxy. If 
the kinematics of the local halo subdwarfs are representative, then the Milky Way's 
dark halo model is oblate with a flattening in the potential of g 0.85, corresponding 
to dark matter density contours with flattening « 0.7. 

The fractional pressure excess for power-law populations is roughly proportional 
to both the ellipticity and the fall-off exponent. Given the same pressure excess, if the 
density profile of one stellar population declines more quickly than that of another, 
then it must be rounder. This implies that the dual halo structure claimed by Carollo 
et al. (2007) for the Galaxy - a flatter inner halo and a rounder outer halo - is 
inconsistent with the virial theorem. 

For the thick disc, we provide formulae for the virial sequences of double- 
exponential discs in logarithmic and Navarro-Frenk- White (NFW) haloes. There are 
good matches to the observational data on the flattening and anisotropy of the thick 
disc if the thin disc is exponential with a short scalelength w 2.6 kpc and local surface 
density of 56 ± 6Mqpc -2 , together with a logarithmic dark halo. Thin discs with long 
scalelengths « 3.5 kpc are disfavoured. Likewise, Navarro-Frenk- White potentials do 
not seem to produce virial sequences matching the thick disc kinematics. 

Key words: galaxies: kinematics and dynamics - dark matter - Local Group Galaxy: 
halo - Galaxy: disc 



1 INTRODUCTION 

The steady-state tensor virial theorem (see e.g. Chan- 
drasekhar 1987) relates the overall geometry of the matter 
distribution to the global kinematics. It is a global result 
that constrains the components of the pressure, kinetic en- 
ergy and potential energy tensors. The potential energy ten- 
sor Wij depends on the flattening, whilst the kinetic energy 
tensor Tij describes the rotational and velocity anisotropy 
properties. 

* E-mail: aagnello@ast.cam.ac.uk, nwe@ast.cam.ac.uk 



Stars and their kinematics provide one of the tools avail- 
able for studies of the morphologies and sizes of galaxies and 
their halos. So, a classic application of the tensor virial theo- 
rem is to the problem of the shape of elliptical and lenticular 
galaxies. Binney (1978, 2005) gave a pictorial representation 
of the plane of ellipticity e versus ratio of rotational to pres- 
sure support v/a for early- type galaxies. He showed that 
if elliptical galaxies have isodensity surfaces that are similar 
ellipsoids, then v/a is completely determined by the velocity 
anisotropy and the flattening. This diagram is an important 
tool in understanding whether the flattening of an elliptical 
galaxy is due to velocity anisotropy or rotation. 
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Although the tensor virial theorem has been widely 
used in galactic astronomy, applications have been largely 
confined to self-gravitating systems - such as the stars in 
bulges or elliptical galaxies — which are moving in a poten- 
tial mainly generated by their own density field. However, 
some populations make a very modest contribution to the 
total potential in which they move. An example is the pop- 
ulation of metal-poor halo stars, which are mere tracers of 
the gravity field of the overwhelmingly dominant dark halo. 
Similarly, the stars of the thick disc move in the combined 
field of the more massive thin disc and dark halo. We shall 
refer to such instances as tracer populations. Here, we pro- 
vide a systematic investigation of the tensor virial theorem 
applied to tracers, including stellar haloes and thick discs. 

Conventional methods of analysis of stellar kinematical 
data on tracer populations are usually based upon phase 
space distribution functions (e.g, Evans, Hafner & de Zeeuw 
1997, Deason, Evans & Belokurov 2011a) or upon the Jeans 
equations (e.g., van der Marel 1991, Battaglia et al. 2005). 
The former is most fundamental as it is equivalent to solv- 
ing the collisionless Boltzmann equation. The method usu- 
ally depends upon an Ansatz as to the likely phase space 
distribution function, with the free parameters fit in a max- 
imum likelihood sense to the kinematical data. The latter is 
less fundamental, as the Jeans equations are moments of the 
collisionless Boltzmann equation. The method depends not 
just on the stellar density and velocity field, but also explic- 
itly on their gradients, which are hard to obtain with noisy 
data and susceptible to deviations caused by contaminants, 
interlopers and substructure. 

The tensor virial theorem is obtained by integration of 
moments of the Jeans equations, and so can be regarded 
as a globally averaged version of them. Although there are 
earlier attempts to use the tensor virial theorem to study 
properties of the Galaxy via the kinematics of tracer stars 
in the literature (White 1989, Sommer-Larsen & Christensen 
1989), they were held up by a serious drawback. Kinematic 
quantities, for example, of halo stars in our Galaxy are only 
known locally, and there is a leap of faith required in ex- 
trapolating from local to global quantities. However, given 
the many large scale surveys of Galactic stars underway or 
planned, it is now worthwhile to examine the tensor virial 
theorem anew for its potentiality. Astrometric satellites such 
as GAIA (Turon et al. 2005), current spectroscopic surveys 
such as SEGUE (Yanny et al. 2009), and forthcoming initia- 
tives, such as ESO-GAIA and Hermes, offer the prospects of 
huge amounts of kinematic information becoming available 
for many stars throughout the Galaxy. Understanding these 
complex and overlapping datasets, and extracting useful in- 
formation from the radial velocities and proper motions, will 
be a major task over the next few years. Global averages 
like the tensor virial theorem may provide a quick way of 
extracting ready and accurate information on the structure 
of the Galactic potential. 

The paper is arranged as follows. In Section 2, we in- 
troduce the tensor virial theorem. Section 3 derives formulae 
for the ratio of equatorial to vertical pressure for scale-free 
tracer populations in spherical haloes. In Section 4, we de- 
velop theorems that extend our results to the cases of flat- 
tened haloes and more general tracer populations. Section 
5 discusses applications to Population II stars in the stel- 
lar halo of the Milky Way. Using observational data on the 



pressure ratio, we develop a new method of constraining 
the flattening of the dark halo. We also show that virial ar- 
guments impose strong constraints on multiple populations 
moving in the same potential and apply our arguments to 
recent claims of duality in the stellar halo. Section 6, dis- 
cusses applications to the thick disc, in particular focusing 
on whether the observed kinematics and flattening are con- 
sistent with motion in a thin disc potential with short or 
long scalelength. Finally, we revisit the stellar halo in Sec- 
tion 7, asking whether it is conceivable that the flattening 
to the potential provided by the thin disc alone reproduce 
the observed flattening and kinematics. 



2 THE TENSOR VIRIAL THEOREM 

Let V be the volume of our survey. For an all-sky survey like 
GAIA, V might be all space. However, this need not always 
be the case. 

Whatever the volume V, the kinetic energy tensor Kij 
is the sum 2Tij + , where 



2 ,3 

pOija. x. 



(1) 
(2) 



Here, p is the stellar density, whilst Vi is the velocity field 
and afj is the velocity dispersion tensor. Here, an overbar 
denotes an average over the velocity distribution so that 



cry = (vi -Vi){vj 



(3) 



We note that Tij is the kinetic energy in ordered, rotational 
motion (often negligible for stellar halos but substantial for 
thick discs) , whilst riy is the pressure tensor. The potential 
energy tensor is 



<9$ , 3 
px l - — d x 

V OXj 



(4) 



where $ is the gravitational potential. 

To obtain the relationship between these tensors, we 
start from the steady-state Jeans equations, which read 



d 
dx k 



■pv k Vj 



dxj 



(5) 



Multiplying by Xi, integrating over the volume V and mak- 
ing use of the divergence theorem gives the steady-state, ten- 
sor virial theorem (e.g., Chandrasekhar 1987, Binney 1981) 



where the surface term is 



Sij — 



XipVjVk.dSk 



(6) 



and dV is the bounding surface. If the integration volume V 
is taken as all space, then the surface term can be discarded 
provided that the pressure pvjvk falls off more quickly than 
r" 3 . 

For an axisymmetric density and potential, we will often 
find it convenient to use cylindrical polars (R,(f>,z). It is 
useful to sum up all the kinetic energy or pressures in the 
equatorial plane and refer to 



K RR . = K xx + K y 



11/ 



U xx + IL 



(8) 
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Analogously, we can also introduce 



Wrr = W xx + Wy 



pRd R $<fx 



and Srr — S X x + S yy . 

The great advantage of virial methods is that integrated 
quantities (such as total kinetic energy of the stars in the 
survey volume) can in general be evaluated robustly even 
with quite sparse data. A price to be paid is that the sur- 
face terms may not necessarily vanish and may have to be 
evaluated from the data. 



3 TRACER POPULATIONS IN SPHERICAL 
DARK HALOES 

For the moment, the gravitational potential is assumed to 
be dominated by dark matter. The stars comprise a tracer 
population. This assumption is fine for stellar haloes which 
move in the dominating dark matter potential. For example, 
the Milky Way Galaxy's stellar halo has a total mass of 
~ 4 x 10 8 M© (e.g., Bell et al. 2008), which is four orders of 
magnitude smaller than the dark matter halo (e.g., Watkins 
et al. 2010). 

3.1 The Pressure Ratio 

Let us take the potential to be spherically symmetric (i.e. 
&(x, y, z) = $(r)) and the stellar density to be axisymmetric 
(i.e., p(x,y,z) = p(R,z)). Let the integration volume V be 
over all space and let the surface terms vanish, so that the 
tensor virial theorem can be manipulated to yield: 



Krr 

K zz 



pR 2 *' 3 



d a x 



(9) 



d*x 



where <£>' = d$/dr. The quantity Krr/K zz is the ratio of 
the equatorial or in-plane kinetic energy to the vertical ki- 
netic energy. 

Stellar haloes are only weakly rotating (e.g., Deason et 
al. 2011a), so it is to fair to assume that there is no mean 
streaming motion (T;j = 0). Then, Krr/ K zz is the same as 
TIrr/TI zz , the ratio of pressure in the equatorial to vertical 
directions. For an axisymmetric figure, it is natural to expect 
this to be intimately related to the flattening. 

3.2 Scale-Free Spheroidal Densities 

To begin with, let us assume that the tracer density is scale- 
free and has the form 



m 



(10) 



where (r, 8) are spherical polars and 7 is the power-law fall- 
off. Formally, the scale-free density is non-integrable either 
at the origin or at infinity and so the virial quantities do not 
exist. Nonetheless, their ratio is well-defined and given by 



Wrr Urr 



f(9) sin 8d8 



u z 



f{8) sin 9 cos 2 6d6 



(11) 



As an illustration, we specialise to the case in which the 
tracer density is stratified on similar concentric spheroids as 



( -ni . 2 -2\-7/ 2 
p = po (R + z q ) 



(12) 

with (R, 4>, z) familiar cylindrical polar coordinates. Here, q 
is the axis ratio of the equidensity contours, which is simply 
related to the ellipticity £ = 1 — q. The virial ratio eq (|ll|l is 



Wrr 

W zz 



11/ 



f- 7T 

Jo 



sin 3 6>(sin 2 6 + q~ 2 cos 2 8) 



-7/2 



dO 



IF 



sin 9 cos #(sin 8 + q cos 



- 7 /2 



dO 

(13) 

which is independent of $ 

Some simple and useful results can be obtained in the 
case of modest flattening, q — > 1. Using the shorthand Q = 
q~ 2 , then Taylor expansion yields 



IF 



IF 



1 1 9 

sin 3 8(1 - -7[Q - 1] cos 2 8)d8 

sin 8 cos 2 6»(1 - i 7 [Q - 1] cos 2 8)d9 



(14) 



We can get an expression for the vertical and horizontal pres- 
sure difference by making the definitions (cr 2 z ) = H zz /M 
and (a xx ) = \Hrr/ M , where M is the total mass of the 
tracer population. The angled brackets therefore denote 
mass-weighted averages of the velocity dispersions. We can 
now recast eq (|14[) as 

» 1+ i w _ 1)+ ii&za w _ 1) . (15) 

Bearing in mind that Q = 1/q 2 = 1/(1 — e) 2 , where e is the 
ellipticity, we see that the fractional pressure excess is 

i! HsP»r«-'>»b- (16) 

In other words, for a mildly flattened stellar halo with a 
given ellipticity and power-law density fall-off residing in a 
spherical dark halo, then 

2 

Fractional pressure excess ss — x Falloff x Ellipticity . (17) 

We note that, when (a 2 x ) — > (a 2 2 ), Q and hence q tend to 1 
(i.e. the stellar halo is spherically stratified), and vice- versa. 
This behaviour in the isotropic pressure limit is general for 
spherical dark matter haloes. 

In the limit of extreme flattening, an asymptotic expan- 
sion of eq (|13[1 yields 



n, 



9 



7iog(2<r 2 ) 



if 7 > 3 

if 7 = 3 (18) 



.i(3-7)9 1 - 7 B(^ i ,i) ifl<7<3 



with B(a;b) = F(a)F(&)/F(a + 6) being the Beta function. 
Stellar halo populations typically have 7 > 3 and so the 
pressure ratio Ilflij/rFj, ~ q~ 2 for high flattening. We will 
see later in Section [6] how the presence of a thin disc alters 
this behaviour. 



3.3 Exact Results 

In fact, the pressure ratio in eq (|13[) is analytic in the case 
that power-law fall-off 7 is an integer. It is helpful to have 
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Figure \T\ shows the behaviour of TL XX /TL ZZ valid for any 
spherical potential. It assumes a power-law profile for the 
stellar density stratified on similar concentric spheroids with 
axis ratio q and with density fall-off m varying between 3 
(lower curve) and 5 (upper curve). Populations in the stellar 
halo typically fall off like r~ 3 or r~ 4 , though some variable 
star populations such as RR Lyrae (e.g., Watkins et al. 2009) 
and Miras (e.g., An et al. 2004) can fall off faster. As a 
general rule, the faster the density fall-off, and the greater 
the flattening, then the larger the pressure excess needed to 
support the stellar distribution. 



Figure 1. Virial sequences for spheroidal stellar halos. The pres- 
sure ratio ^nfl/j/IIjjz is shown as a function of the flattening q 
for stellar densities with fall-off 7 = 3 (blue), 3.5 (green), 4 (yel- 
low), 4.5 (red) and 5 (purple). The pre-factor of i means that the 
vertical axis is Ii X x/H zz and ensures that the ratio tends to unity 
as q — > 1 . The analytic form of these curves is given in eqs 1201 1 , 
(EQ and J22} for 7 = 3, 4 and 5. 



the exact expressions. For 7 = 2, we find: 



11/ 



» arctan y/Q-1 - (Q - 1) 3 



n. 



(Q — 1) 2 — arctan VQ — 1 



(19) 



For the case 7 = 3, we have: 

Ilftfl _ (Q- 1)1 - ^arcsinh^Q - "! + (Q- 1)4 



IF 



V / Qarcsinhy / Q _ 



(0-1)- 



(20) 

This formula is equivalent to one derived by White (1989) 
for populations orbiting in an isothermal sphere. 

The 7 = 4 case has already been discussed by Helmi 
(2008). We find: 



n RR _ Q(Q -2) arctan VQ~~1 + (Q - 1) 



n. 



Q arctan y/Q — 1 — (Q — 1) 2 



(21) 



Like us, Helmi assumes that the stellar halo is flattened and 
follows a scale- free power-law density, whilst the dark halo is 
spherical. In fact, Helmi specialises to the particular case of 
a logarithmic potential generating an asymptotically flat ro- 
tation curve. To make progress, Helmi makes an additional 
assumption, namely that the ratio of horizontal to vertical 
velocity dispersions o^/cr 2 does not depend on location. Un- 
fortunately, the final formula - equation (8) of Helmi (2008) 
- contains two mistakes. The first is computational, namely 
the integrals are incorrect. This can be seen by noting that, 
in the limit of a spherical stellar halo (Q — > 1), the veloc- 
ity dispersion tensor must become isotropic within Helmi's 
set of assumptions, a fact that is contradicted by her for- 
mula. The second is more subtle and conceptual. The ratio 
cr 2 /cr 2 between velocity dispersions in the x and z coordi- 
nates cannot be considered uniform, once the density and 
potential profiles are already assigned: this additional as- 
sumption over-determines the problem. 

For the 7 = 5 case, the result is exceptionally simple, 
namely 



n, 



IF 



2Q 



(22) 



The result is exact for arbitrary flattening. 



3.4 Alignment of the Velocity Ellipsoid 

To make practical use of the virial formulae, we must find 
the connection between the pressure ratio and the velocity 
dispersions. This is straightforward if the velocity dispersion 
tensor is diagonal in cylindrical polar coordinates, namely 



IF 



IF 



p{°RR + 0-^)d 3 X 



2 ,3 

pa zz d x 



(23) 



If the velocity dispersion tensor is aligned in spheri- 
cal polar coordinates - as is the case for halo stars (e.g., 
Smith et al. 2009b) - then matters are more complicated. 
The diagonal components of the velocity dispersion tensor 
in cylindrical polar coordinates are now given by 



2 2 ■ 2 2 2 

a RR ~ a rr sm *> + a 80 cos t 

0%<f, = CTgg = (1 — P) a1 r 
2 2 2 r\ 1 2 ■ 2 / 

a zz = a rr cos a + Ogg sm 6 



/3 cos 9) a rr 



= (1 - /3sin 
where j3 is the anisotropy parameter, 



This means that 



11/ 



IF 



pa^ r (2-/3 + /?cos 2 6»)d 3 a; 



p<7Vr-(l — sin 2 0)d' i x 



(24a) 
(24b) 
(24c) 



(25) 



(26) 



If P is constant, this simplifies to 

Urr _ (2~P)(a 2 rr ) -P{al r 



n Z2 



(27) 



where angled brackets denote mass- weighted averages. If the 
tracer flattening is known and the variation of the veloc- 
ity dispersions is estimated over the survey volume, eq (|27[1 
yields virial sequences that do not involve the potential $. 
Thus, we can constrain the potential by intersecting the 
Wrr/Wzz sequences with the pressure ratios. Nonetheless, 
practical application of eq (|27|) is limited by the fact that 
the global variation of the dispersions with position is not 
well understood, though this will be remedied after the next 
generation of large-scale surveys of the Milky Way (such as 
GAIA). 
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Logarithmic potential NFW potential 



R 




Figure 2. Virial sequences for the stellar halo density law found by Deason et al. (2011b) and given in eq. (|3Q l> - The pressure ratio 
2 ^rr/TIzz is shown as a function of flattening ratio q/g, where q is the axis ratio of the stellar density contours and g is the axis ratio 
of the dark matter equipotentials. The dark matter potential has a logarithmic form generating a flat rotation curve (left) and a NFW 
dark matter halo (right). The coloured lines refer to different values of the dark matter scalelength from 5 kpc (lower red curve) through 
15 kpc, 35 kpc, 55 kpc, 75 kpc to 95 kpc (upper blue curve) and taking the break radius fixed at R\, = 27 kpc (Deason et al. 2011b). 
The green and blue bands give the range in pressure ratios found by global extrapolation of the local results of Kepley et al. (2007) and 
by Smith et al. (2009a). 



4 TWO THEOREMS 

4.1 The Flattening Theorem 

If our formulae were restricted to spherical dark haloes, they 
would have limited applicability. Here, we derive a remark- 
able theorem that extends our results to the case of flattened 
tracer populations in flattened dark haloes. 

Let the tracer density and total potential satisfy the 
following requirements 



p(R,qz) = p(qz,R), 



<f>(R,gz) = $(gz,R), (28) 



for some q and g. We say in this case that p and $ are 
stratified on constant flattening surfaces. The above defini- 
tion is equivalent to requiring that p (or $) be a symmetric 
function of R and z/q (or z/g). The most commonly used 
examples are spheroidal distributions, such as 



$ = < l > (m), 



m 2 = R 2 + z 2 g 



(29) 



There are a number of popular halo models of this form, 
including Binney's (1978) logarithmic potential and the 
power-law models (Evans 1993, 1994). The dark matter den- 
sity may be cored or cusped. The flattening, or the minor 
to major axis ratio, of the dark matter equipotentials g is 
rounder than the dark matter isodensity contours. Similarly, 
the spheroidal power-law densities with fall-off 7 and axis 
ratio q studied in the previous section belong to this class, 
though there are many other possibilities for the tracer den- 
sity as well. 

In what follows, it is convenient to specialize to 
spheroidally stratified potentials with flattening g. The ten- 
sor virial theorem becomes: 



IT 



IL 



where $' = d$/dm. 



pRdR$(Tx 
pzd z <&A 3 x 



pR 2 $' 3 



d°x 



pg- 2 z 2 &, 3 



(30) 



Making the substitution 2 = z'g, writing r' 2 = x 2 + y 2 + 
z' 2 and then dropping the primes, the integrals on the right- 
hand side of eq (|3Up are converted into the exact form as 
those on the right-hand side of eq ©, with the only change 



p(R,z/q) -> p{R,gz/q). 



(31) 



This means that all our preceding results for spherical haloes 
hold good in the flattened case, but with the replacement 



q -> q/g- 



(32) 



This considerably extends the scope of our earlier results. 

More formally, our theorem can be stated as follows. // 
the dark matter potential $ is homoeoidally stratified with 
flattening g and the tracer density p has a constant flatten- 
ing q, then the virial ratio Wrr/W zz depends on the two 
flattenmgs just through their ratio q/g. 

As a corollary, for scale-free densities with constant flat- 
tening the virial ratios depend on q/g and the power-law 
exponent, but not on the radial form of the potential, as 
already noticed for spherical potentials. This can be seen 
by using spherical polar coordinates and noting that the ra- 
dial integrals (containing $) factorize and cancel top and 
bottom. The angular integrals do depend on the density 
stratification and the behaviour of Wrr/Wzz, and thence 
TIrr/TIzz, is controlled by those: 



Uzz 



p(sin0, (g/q) cos 6) sin 3 060 



p(sin 9, (g/q) cos 9) cos 9sm9d9 



(33) 



As a consequence, once a virial ratio has been computed 
for a spherical dark matter potential, the flattened analogue 
follows in a straightforward manner. Providing the poten- 
tial remains spheroidally stratified, any halo model such as 
isothermal or Navarro-Frenk- White (1998) or Sersic (1963) 
is fine. This is an extension of a famous and elegant result 
from potential theory due to Roberts (1962), which states 
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that ratios such as Wrr/W zz depends only on the axial 
ratios for homoeoidal distributions, and not on the radial 
profiles. 

The ratio Wrr/W zz does depend on the form of 4>, if 
more complex profiles are prescribed, for example, account- 
ing for the contribution of the stellar disc, or the differential 
flattening of the dark matter halo and the self-gravity of 
the stellar halo. However, although analytic progress is dif- 
ficult for complex profiles, the relevant integrals are easily 
evaluated numerically. 



4.2 The Conjugateness Theorem 

As shown above, there is a class of test densities (like scale- 
free ones) yielding virial ratios that are independent of the 
radial behaviour of $. Similarly, we can ask whether there 
are potentials giving virial ratios that are independent of the 
radial stratification of p. If a density and a potential belong 
to these classes and give the same virial ratios, we refer 
to them as conjugate. We shall demonstrate that scale-free 
potentials are conjugate to scale-free tracer densities. 

Consider a scale-free power-law potential describing the 
dark matter 



$(r) = < 



2 



if a 0, 



(34) 



vq log r if a = 



where a is a constant. Now, the circular velocity in the 
equatorial plane varies like r~ a ^ 2 . So models with a < 
have falling rotation curves, models with a > have ris- 
ing rotation curves. When a = 0, the model is the familiar 
isothermal sphere with a perfectly flat rotation curve vq. 
Now suppose that the tracer density has the form 



(35) 



Note that this density profile is not necessarily scale-free, 
but could be of exponential or Gaussian form. Then, the 
ratio of equatorial to vertical pressure is 



11/ 



[h{9)f +a sin 3 6A6 



n. 



(36) 



[h{6)f +a sinO cos 2 6d6 
But, this is the same virial ratio as would be given by 

$ = $(r), p=^- (37) 



with 7 = a + 3. The two combinations are conjugate. 

The conjugateness theorem then reads: A scale-free 
tracer density with fall-off a + 3 in an arbitrary spherical 
potential is conjugate to a tracer density in a scale-free po- 
tential with fall-off a. To give a practical application of this 
result, we might want to know the pressure ratio of a tracer 
population with 



po exp 



2 . 2 
+ z q 



-2\l/2 



/Rt 



(38) 



in a galaxy with a flat rotation curve (a = 0). The density 
law might describe a thick disc with horizontal scalelength 
-Rthick and vertical scaleheight q -Rthick. By the conjugateness 



Datasct 


Dark Halo 


n RR /u zz 


l/9 


9 




Model 




K07 


NFW 


5.28 ± 0.52 


0.52 ±0.06 


1.2 ± 0.14 


K07 


Logarithmic 


5.28 ± 0.52 


0.55 ±0.06 


1.1 ± 0.12 


S09 


NFW 


3.26 ± 0.10 


0.70 ±0.03 


0.85 ±0.04 


S09 


Logarithmic 


3.26 ± 0.10 


0.73 ±0.03 


0.82 ±0.04 



Table 1. Results for the ratio of stellar halo to dark halo flat- 
tcnings q/g for different assumed halo models using the Kepley 
et al. (2007) and the Smith et al. (2009a) datasets. The value of 
the dark halo axis ratio g is given assuming that the stellar halo 
flattening is q = 0.6 as found by Deason et al. (2011b). 



theorem, this is just the same as the pressure ratio for a 
spheroidal population with 7 = 3, which is given in eq (1201) . 

We finally remark that, although the potential and the 
density may be chosen as scale-free, they cannot both be 
so chosen as the virial ratios are then undefined. White's 
(1989) results for spheroidal tracer populations orbiting in 
the singular isothermal sphere hold good providing the core 
radius of the tracers R c does not vanish. Sommer-Larsen 
& Christensen (1989) give formulae for the virial ratios of 
scale-free populations in scale-free potentials, in which the 
core radius of the dark matter Rdm is also zero. However, 
the virial ratios depend on the manner in which the limit 
R c — > and -Rdm — > is approached, and so their results 
must be treated with caution. 

More precisely, given a power-law density and potential 
with core-radii R c and Rdm respectively, the virial ratios can 
be easily shown to depend just on the ratio R c /Rd m - The 
limits R c — > and Rdm can be approached by keeping 
the ratio fixed. Since different R c /-Rdm give different virial 
ratios, the final limit depends on the path by which the limit 
is approached, and so is not defined. 



5 STELLAR HALOES 

5.1 The Flattening of the Milky Way Dark Halo 

The flattening in the stellar density q and the pressure ratio 
T1rr/TI zz are in principle accessible directly from starcounts 
and kinematical measurements of halo stars. The flattening 
in the dark halo density g is not. The Flattening Theorem 
therefore provides a completely new method to study the 
dark halo flattening. 

Deason et al. (2011b) traced the stellar halo of the Milky 
Way using blue stragglers and blue horizontal branch stars 
as tracers. They found that the halo has a broken density 
profile given by 



p(R,z) = < 



Po 



V-R 2 + z 2 /<3 2 
^/R 2 + Z 2 /q 2 



z 2 /qi < R b , 



z 2 /q 2 R b . 



(39) 

This law has an explicit break radius Rb, at which the power- 
law profile changes from a shallow inner slope of 2.35 to a 
steeper outer one of 4.6. Deason et al. (2011b) presented 
evidence that Rh = 27 kpc and that q = 0.6. 
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Fig. [2] shows curves of pressure ratio for the stellar den- 
sity law (|39[) as a function of q/g, where q is the axis ratio 
of the stellar density contours, whilst g is the axis ratio of 
the dark matter equipotentials. Now, the density law is no 
longer a homogeneous function of coordinates, and so our 
results do depend on the radial profile of the gravitational 
potential. The left panel gives results for a logarithmic po- 
tential 

2 

$ = ^log(i? dm 2 +m 2 ), m 2 = R 2 + z 2 g~ 2 . (40) 

The right panel shows the outcome for the same stellar 
density, but now in a spheroidally stratified Navarro- Frenk- 
White potential 

47rp n f w i?dm 3 . / 1 , m \ , . 

$ = ^ log 1 + , (41) 

m V "dm / 

which has an everywhere positive density providing g > 
1/V2. The different coloured lines show the effect of varying 
Rdm between 5 kpc (lower curve) and 95 kpc (upper curve) , 
using a fixed break radius Rb = 27 kpc as advocated by 
Deason et al. (2011b). 

The stellar kinematics of halo stars in the local neigh- 
bourhood have been measured many times. Specifically, we 
assume that the locally calculated ratio of stellar velocity 
dispersions is a good approximation to the density-weighted 
globally averaged ratio. Even though the velocity dispersions 
of course change with position, the local ratio can still be 
an excellent approximation to the global one provided the 
anisotropy does not vary dramatically with position. No- 
tice that we are not explicitly assuming either spherical or 
cylindrical alignment (which in any case coincide locally), 
although we are arguably making a stronger assumption in 
taking local ratios as representative of global ones. Of course, 
this limitation to our calculation will vanish once surveys 
such as GAIA have been completed and provide a global 
picture of how the velocity dispersion varies for components 
in the Galaxy. 

Kepley et al. (2007) assembled a sample of 221 halo 
stars within 2.5 kpc of the Sun with distances, radial ve- 
locities and proper motions. They found (ctrr, cr^, cr zz ) = 
(175±8, 110±6, 84±4) kms -1 around means of {vr,v^,v z ) — 
(-4 ± 11, -23 ± 8, -1 ± 6) kms" 1 . This gives Urr/U zz = 
5.28±0.52, which is represented as the green horizontal band 
in Fig. [2] As the Kepley et al. data are drawn from a va- 
riety of disparate sources, the errors are appreciable. Smith 
et al. (2009a) used a sample of ~ 1700 halo subdwarfs ex- 
tracted from the light-motion curve catalogue of Bramich et 
al. (2008) of Sloan Digital Sky Survey (SDSS) Stripe 82. The 
halo subdwarfs were extracted using a reduced proper mo- 
tion diagram. This is a larger and more homogeneous sam- 
ple of halo stars centered on R = 8.70 kpc and z — —2.41 
kpc, and so below the Galactic plane. Smith et al. found 
{o RR ,a^,a zz ) = (138 ± 2, 82 ± 1,89 ± 1) kms" 1 around 
means of (v R ,v,p,v z ) = (9 ± 3, 2 ± 2, -1 ± 2) kms" 1 . This 
gives Hrr/H zz = 3.26 ±0.10, which is shown as a light blue 
horizontal band on Fig. [2] 

Table [l] summarizes the results for the different datasets 
and halo models. In particular, the Kepley et al. (2007) 
dataset is consistent with a spherical dark matter halo. The 
Smith et al. (2009a) dataset is consistent with oblate dark 
halo models with a flattening in the potential of g ~ 0.85, 





Short Disc 


Long Disc 


6 


[0.92,2.25] 


[0.25,0.56] 


d 


[0.49,0.81] 


[0.78, 0.95] 




[0.3,0.7] 


[0.7,1.4] 



Table 2. Parameters for the the disc-halo decompositions 
adopted in Section \E\ assuming that the dark matter potential 
has a logarithmic profile. 



corresponding to a flattening in the dark matter density con- 
tours of ~ 0.7. Both results assume that there is no con- 
tribution from the thin disc potential and so this analysis 
provides lower limits to the axis ratio q. We shall revisit the 
effects of the thin disc in Section [7] 

5.2 The Dual Halo 

Even though the result given in words in eq (|17|l is a simple 
one, it can still be made to do some work. Here, we shall 
show that it implies that the dual halo structure advocated 
by Carollo et al. (2007, 2010) is unphysical. 

Carollo et al. (2007) first argued for a two component 
halo on the basis of different density profiles, metallicities 
and stellar orbits. They claimed evidence for the existence 
of a metal-rich, inner halo ([Fe/H] ~ -1.6) and a metal-poor 
outer halo ([Fe/H] ~ -2.2) The inner halo is supposed to 
have a modest net prograde rotation, together with an oblate 
shape with an axis ratio of q ~ 0.6. The outer halo has net 
retrograde rotation with a much more spherical distribution 
with an axis ratio of q ~ 0.9 to 1.0. 

Carollo et al. (2010) supplied further structural and 
kinematical analysis of the dual halo. The inner halo has 
a density profile that falls off like p ~ r -3- 17 ±0-20 T ^ ye _ 
locity dispersions (ctrr, 0$$, a zz ) = (150 ± 2, 95 ± 2, 85 ± 2) 
kms - around means of (vr, v^,v z ) = (3 ± 2, 7 ± 4, 3 ± 1) 
kms -1 . This gives Krr/K zz = 4.40. for the inner halo. 

The outer halo has a density profile that falls off like 
p ~ r - 1 -79±o.29 jj. nag ve io C ity dispersions (grr, cr<f,<t>, a zz ) = 
(159 ± 4, 165 ± 9, 116 ± 3) kms -1 around means of 
(vr,v^,,v z ) = (-9 ± 6,-80 ± 13, 2 ± 4) kms -1 . This gives 
Krr/K zz = 4.40, for the outer halo, surprisingly close to 
the value obtained for the inner halo. 

Given that the fractional pressure support is the same 
for the inner and outer halo, we see that the claimed Batten- 
ings and power-law fall-offs are inconsistent with the virial 
theorem. If the density profile of the inner halo declines more 
quickly than that of the outer halo, then it cannot be flatter 
as well. Put simply, the claimed shapes, profiles and kine- 
matics of the inner and outer halo are inconsistent with the 
structures moving in the same steady-state potential. If the 
outer halo is rounder than the inner halo, its density profile 
must fall off more quickly than the inner one. 

It is useful to summarise the assumptions in this con- 
clusion. The result does assume that q, the flattening of the 
stellar halo, is constant. There is good evidence for this as 
investigations of changes of ellipticity with radius have typi- 
cally found none (e.g., Sluis & Arnold 1998; Sesar et al. 2011; 
Deason et al. 2011b). The result also assumes that g, the flat- 
tening of the dark matter equipotentials, is constant. Whilst 
there is no direct evidence for this, it certainly seems reason- 
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1.41 
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1 9 


3 


8 


1 57 


10 


1.7 


2.5 
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1.37 


10 


6.2 


3.0 


•3.6 


1.55 


12 


.3.7 


2.5 


8.1 


1.37 


12 


3.8 


2.5 


7.5 


1.42 


12 


1.7 


3.0 


2.8 


1.32 


10 


0.9 


2.5 


3.0 


1.42 


10 


1.5 


3.0 


1.1 


1.39 


40 


2.0 


1.5 


0.6 


1.50 


200 


0.1 


2.1 


1.10 


1.46 


200 


0.1 


2.8 


0.28 


1.29 


200 


0.2 


3.1 


0.20 


1.45 


200 


0.2 


3.5 


0.15 


1.50 



Table 3. Sample parameter values for the decompositions with 
an exponential disc with scalelength iJthin and NFW potential 
with scalelength i?dm an d concentration c . The rotation curves 
are shown in Fig. \3\ All the rotation curve fits take the surface 
density of the thin disc as 56 Mq pc — 2 at the solar radius of 8.0 
kpc. 



350 - 



300 - 




5 10 15 20 

R (kpc) 



Figure 3. The rotation curves of model galaxies with NFW 
haloes and exponential thin discs with the parameters listed in 
Table [3] The data arc taken from Sofuc ct al. (2009). The fit is 
performed using the datapoints beyond 3 kpc, for which the con- 
tribution of non-circular motions caused by the bar is insignifi- 
cant. 



able to investigate such simple models first of all. Although 
we have used the asymptotic formula (|17[) in our argument, 
the problem for Carollo et al.'s model becomes worse if the 
full result is used. As can be seen from the curves in Fig. [TJ 
a greater pressure excess is needed at larger flattenings than 
implied by a linear scaling with ellipticity. 



6 THICK DISCS 

Our methods are also applicable to thick discs, and can give 
some simple results. As thick discs are predominantly ro- 
tationally supported, it makes sense to consider the total 
kinetic energy rather than the pressure contribution alone. 
Mass ordering arguments suggest that we discard the con- 
tribution of the bulge in the virial relations, even if it is 
important in shaping the innermost parts of the rotation 
curve. 



Long Disc 

2J — i — 1 — 1 — 1 — i — 1 — 1 — 1 — i — 1 — ' 



2.0 - 




-1.2 -1.0 -0.8 -0.6 -0.4 

Log 10 (?) 

Figure 4. Virial sequences for the thick disc. The equatorial to 
vertical kinetic energy ratio is shown as a function of the flattening 
q. The thick disc stars move in the potential of a long thin disc and 
logarithmic halo with relative mass ratio fj, = 0.7. Each diagonal 
stripe corresponds to fixing b = i?thick/^dm an d varying d = 
Kthin/^?thick as in Table [2] The vertical and horizontal bands 
indicate the obscrvationally preferred values of the flattening and 
the kinetic energy ratio. Notice that all the models lie well below 
the observational window. [Blue curves are for b = 0.2, purple 
for b = 0.25, green for b = 0.3 and orange for b = 0.35. If fi is 
increased (decreased), the virial sequences move down (up) in the 
figure.] 



6.1 Formalism 

The density profiles of galactic discs tend to be exponen- 
tial or Gaussian - more complicated than the homogeneous 
functions used to represent stellar haloes. We shall assume 
that the thick disc has a double-exponential form: 

^ oexp (-^) exp ( _ 5) ^ 

and so the vertical scaleheight — gJ?thick, where q is as 
before the axis ratio of the density contours. 

Thick discs move in the potential of both the dark halo 
and thin disc. The latter can be represented by a razor- 
thirQ exponential disc with scale- length i?thm, and thus with 
potential: 

2 f°° J (kR)e- k ^ AU 
$ disc (i?, z) = E ^ [1 + Wfc2]3/2 dfc. 

(43) 

The ratio of the kinetic energy tensor components takes the 
form 

/ pRd R ($ diBC + $hE,io)d 3 a; 
^ = ~f ' 

and so depends on both the thin disc and halo potentials. 
By inserting the thick disc density (|42[) and the thin disc po- 
tential (|43[1 into (|44|l , we can find that the disc contributions 



1 We have calculated the correction to our virial terms by dis- 
carding the assumption of zero-thickness of the thin disc. We find 
it to be always < 7 % if the thin disk scaleheight is < 350 pc for 
realistic values of the local surface density (0.3 < fi < 0.7). 
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are 

Dr oc 



j pRd R $di a cd 3 x 



3k 2 



(l + gfe)(l + fe 2 ) 5 /2(l + fc2 d 2)3/2 



D z oc / pz9,<3?di sc d x 



k 



(1 + gfc) 2 (l + fc 2 ) 3 / 2 (l + fc 2 d 2 ) 3 A 



dfc 



■dk 



(45) 



(46) 



where d = -Rthm/^thick- The choice of a double-exponential 
disc enables the quadruple integrals to be reduced to a single 
one using formulae [6.623] of Gradshteyn & Ryzhik (1965). 
This simplification is not possible for general thick-disc den- 
sity profiles. 

Now, let us suppose that the halo is the spherical limit 
of the logarithmic model 

$o log (r 2 + R dm 2 ) . 



$halo 

Then the halo contributions are 

.X 3 exp-(x + £) 



H R = b 

H z = bq 



dX /o dC ?+^W 
, c x£ 2 exp-(x + £) 



(47) 

(48) 
(49) 



where b — -Rthick/-Rdm- With some careful algebra, we can 
recast eq. (|44"|) as 



K RR 1 2H R + nD R 



K zz 

where the constant p is 



M = 



2H Z + pD z 



27rE -RthiE 



(50) 



(51) 



$0-Rdm 

In this case, the flattening is not fixed and there is no sim- 
ple dependence as was the case in the previous sections. The 
flattening in the total potential is controlled by the p param- 
eter, which mediates the contributions of the razor-thin disc 
and the round dark halo. 

As an alternative model for the dark matter potential, 
we also examine the spherical NFW profile, for which the 
gravitational force component is 

47rGp n f w i?dm 3 



log 1 + 



Rdn 



r+R&n 



ATlGpnlwRdn 



(52) 



The disc contributions remain the same, but the halo ones 
become 

X 3 exp-(x + 



H R = 
H z = 



, dx / dC (x 2 +<? 2 e 2 ) 3/2 



HWx 2 we) ,(53) 



t xC 2 exp-(x+C) a 



and the constant p is now 



A 4 = 



So^thin 2 



Hbyj x 2 we)m) 



(55) 



pnfw Rd 

In the limit of a massless thin disc p — > 0, the ratio 
Krr/K zz is roughly oc q~ 2 in the highly flattened limit. 
This can be deduced from eq. (|18[) using the conjugateness 



theorem. Letting the disc become massive, or equivalently 
introducing a non-vanishing p, then the rise in pressure 
deficit for small q is shallower, approximately oc q^ 1 . This 
is understandable because the presence of a thin disc flat- 
tens the total potential, and thus a smaller pressure deficit 
is required to attain a given tracer flattening when a disc is 
present. 



6.2 Parameter Values 

To make progress, we will need to choose astrophysically 
motivated values for our thick disc models. Here, we are 
limited by the fact that the thick disc of our Galaxy has not 
been fully mapped out. We shall see that with the existing 
data, a rather wide range of parameters is plausible a priori. 

De Jong et al. (2010) used SEGUE photometry to find 
the scalelength of the thick disc i?thick = 4.1 ± 0.4 kpc and 
the scaleheight Zh = 0.75±0.07 kpc. This gives the flattening 
or axis ratio as 



Q = 



0.75 ±0.07 
4.1 ±0.4 



0.18 x (1±0.19) 



(56) 



T his result is supported b y the s urvey of 34 edge-on galaxies 
bv lYoachim fc DalcantorJ (|2006l ) who found that thick discs 
typically have q = 0.18 ± 0.04. 

The velocity dispersions in the thick disc vary 
with radius, perhaps in a roughly exponential man- 
ner (|Lewis fc Freeman! Il989h . The radial variation has not 
been observationally determined in full, but this does not 
harm us here. What is important in virial arguments is the 
ratio of the global averages of horizontal to vertical kinetic 
energy. This ratio can be well constrained even though the 
variation of velocity dispersion in the thick disc is not yet 
fully mapped out. In particular, if the velocity dispersions 
are proportional to ^fp and therefore vary exponentially, 
then the ratio of the local averages of horizontal to vertical 
kinetic energy is equal to the global averages. This provides 
some security that using comparatively local measurements 
of the thick disc velocity dispersions leads to meaningful ki- 
netic energy ratios. 

Carollo et al. (2010) explored a local volume within 4 
kpc of the Sun and reported that the velocity ellipsoid of 
the thick disc is (ctrr, a^, a zz ) = (53 ± 2, 51 ± 2,35 ± 1) 
kms -1 around means of (vr, v^jVz) = (3 ± 2, 182 ± 2, ± 1) 
kms -1 . This gives K RR /K ZZ = 31.5 x (1±0.06). There are 
also v alues available from the work of lCasetti-Dinescu et al.l 
|201ll ). who used red clump stars between 5 and 10 kpc 
of the Galactic Center and between 1 and 3 kpc from the 
Galactic plane. Stars closer to the Galactic plane are not 
used to minimize contamination from the thin disc. For the 
thick disc, they found (<jrr, u^, o zz ) — (70±4, 48±8, 36±4) 
kms -1 around a rotation velocity of = 175 ± 2 kms -1 . 
This yields K RR /K ZZ = 29.2 x (1 ± 0.23). Although the 
values of the dispersions found in the two investigations are 
not the same, the pressure ratio is much more robust, which 



2 There is some t entativ e evidence in Figure 12 of 
Casctti-Dincsc u et al. (2011) that the vertical velocity dis- 
persion of the thick disc does follow a different law as compared 
to the horizontal velocity dispersions. Nonetheless, the error bars 
on these measurements remain quite large (see e.g. Table 3 of 
their paper). 
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Short Disc Short Disc and Bulge 
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0.5 - - 0.5 - 



-1.2 -1.0 -0.8 -0.6 -0.4 -1. 2 -I.0 -0.8 -0.6 -0.4 

L°g 10 (ijO Log 10 (<?) 

Figure 5. Virial sequences for the thick disc. Left: The thick disc resides in the potential of a short thin disc and logarithmic halo with 
relative mass ratio /i = 0.7. Each diagonal stripe corresponds to fixing b = i?thick/^dm ar >d varying d = fithin /-^thick as m Table[2] Now 
all the models cross the observational window, which is the intersection of the vertical orange and horizontal green stripes. Right: The 
thick disc stars now move in the potential of a short thin disc, a Sersic bulge and a logarithmic halo. Notice that the effects of including 
the bulge arc small in the flattening regime indicated by the observationally preferred orange vertical band. [Blue curves are for b = 1.35, 
purple for b = 1.63, green for b = 1.91 and orange for b = 2.25. If fi is increased (decreased), the virial sequences move down (up) in the 
figure.] 



adds weight to our argument that the global ratio is well 
constrained. 

The halo and thin disc parameters are of course strongly 
correlated, as the combined force field must generate the 
Galactic rotation curve. Sofue et al. (2009, 2010) show that a 
decomposition into a de Vaucouleurs bulge, exponential disc 
with spiral arms, and isothermal dark halo gives values of 
the halo core radius Rdm ~ 8— 15 kpc. Sofue et al.'s thin disc 
has a local surface density of 87M0pc~ 2 and a scalelength 
■Rthin = 3.5 kpc - some what larger than is found by most 
other investigators (e.g.,lFux fc Martinet! 1 1994 lOihal l200ll : 
ICasetti-Dinescu et al.ll201l[ ) - and accordingly we refer to it 
as the long disc model. 

On the other hand, iHolmberg fc Flvnnl l|2004 ) used Hip- 
parcos data to measure the local disc surface density as 
56 ± 6Mapc~ 2 , lower than t h at est imated by ISofue et al.l 
(|2010l ). ICasetti-Dinescu et all (|201ll ) used kinematical and 
photometric data to infer the thin disc scale length as 
iithin = 2.6 ± 0.4 kpc. These observations suggest that the 
thi n disc may b e less massive and smaller than advocated 
by ISofue et al.l (|2009h , and accordingly we refer to this as 
the short disc model. To fit the rotation curve, this latter 
model requires lower values of -Rdm, accordingly in the range 
~ 1.5 — 3 kpc. Of course, this is an underestimate of the true 
scalelength of the dark matter, as we have not explicitly in- 
cluded a bulge. With these two choices (short versus long 
disc), the ratios b, d and fi are allowed to vary in the ranges 
summarized in Table [2] Under the assumption of a logarith- 
mic halo, the mass ratio /i in the range ~ 0.7—1.4 (long disc) 
and ~ 0.3 — 0.7 (short disc) to provide a reasonable match to 
the rotation curve. Similarly, the ranges of b = -Rthick/Rdm 
are constrained to lie in the range 0.25 < b < 0.56 (long) and 
0-92 < b < 2.25 (short). From the scalelengths of the thin 
and thick discs, we can deduce the range of d — -Rthm/Rthick 
as 0.78 < d < 0.95 (long) and 0.49 < d < 0.81 (short). 

Matters are more complicated if a NFW profile is cho- 
sen for the da rk matter potentia l. Here, we fit the rotation 
curve data of ISofue et all (|2009h with an exponential disc 
and NFW halo, as shown in Fig. [3] Even keeping the local 



NFW potential 




Figure 6. As Fig. \E\ but for a NFW rather than the logarith- 
mic halo. Plotted are all the virial sequences for the models with 
rotation curves displayed in Fig. [3] Different colours stand for dif- 
ferent values of b = i? t hick/-Rthin, but the sequences are so near to 
one another that listing the colour correspondence is superfluous. 



disc density as 56Mqpc -2 and i? t hin > 2.4 kpc, there is quite 
a large range in which the other parameters may vary. We 
list the parameters of the models in Table [3] although we 
shall show that their virial sequences are not too different. 



6.3 Short versus Long Thin Discs 

Our method is to take i? t hick = 4.1 ± 0.4 kpc from de Jong 
et al. (2010) and -RtMn either short at 2.6 ±0.4 kpc (Casetti- 
Dinescu et al (2011) or long at 3.5 kpc from Sofue et al 
(2009). Then a fit to Sofue's data - as judged by the like- 
lihood falling to half its maximum value - is used to find 
plausible ranges of d and /j, . Virial sequences are found by 
fixing 6, and fi = 0.7, then varying q to get a single line 
across the plane, and finally d (within the range found from 
the rotation curve) to get a band of finite width rather than 
a single line. 

The panels of Figs [4] and [5] show such virial sequences. 
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They give the thick disc kinetic energy ratio as a function 
of flattening q. The allowed observational range is shown by 
vertical and horizontal stripes, so that the intersection gives 
the region compatible with the data. Specifically, the verti- 
cal orange strip e corresponds t o the flattening of the thick 
disc inferred bv lde Jong et all (|201fJ ). whilst the horizontal 
green stripes corres pond to the kinetic energy ratios from 
ICarollo et ail (|201Ch and ICasetti-Dinescu et all |201ll ). 

Plotted on Figs 0] and [5] are the virial sequences of 
the long and short discs respectively. We see that the short 
disc does indeed give a better match to the observations 
than the long disc, even when the fi parameter is taken to- 
wards the high end of its possible range. By contrast, the 
virial sequences for the long disc yield pressure ratios that 
lie appreciably below the estimated ones, for all the relevant 
values of /j, parameter. 

The effects of incorporating a bulge are shown in the 
right panel of Fig [S] The fit to the rotation curve is per- 
formed again, using a short thin disc, an isothermal dark 
halo and a spherical Sersic m = 3 bulge, given by the de- 
projection of 

£(i?) = E exp (-K(R/Re)i^ (57) 

Here, £ is the projected bulge density on the sky, R is pro- 
jected radius and R c is the effective radius and k is a con- 
stant. The values of the parameters are chosen to give a total 
bulge mass of ~ 10 10 M Q and a maximum contribution to 
the rotation curve of ~ 100 kms -1 . The parameter values 
are suggestive rather than definitive, as we merely wish to 
ascertain the importance of the effect. In fact, we notice that 
the virial sequences do not change much in the observational 
regime - the left and right panels of Fig. [5] are very similar 
for q > 0.1. This is understandable, as the models with or 
without a bulge reproduce the same overall rotation curve, 
albeit using different components. The virial quantities do 
not change much, as they depend on the total potential. 

We conclude that the flattening and kinematics of 
the thick disc are consistent with stellar motion in a po- 
tential generated by a thin disc with a short scalelength 
~ 2.6 kpc, togethe r with the normalisation suggested by 
iHolmberg fc Flvnnl (|2004 ). In particular, thin discs with 
longer scalelengths ~ 3.5 kpc are necessarily more massive if 
they are to reproduce the local surface density. The stronger 
gravity field towards the equatorial plane then causes the 
thick disc to be flatter than is measured, given the observed 
kinematics. 

We briefly show the differences made by using a NFW 
halo rather than a logarithmic one in Fig. [5] Here, all the 
virial sequences corresponding to the parameters in Table [3] 
lie below the observational estimates for flattening versus 
pressure excess. The reason for this is as follows. For cusped 
potentials, the virial terms are dominated by contributions 
from near the central cusp. Near the very centre, the com- 
bined equipotentials due to the thin disc and spherical halo 
are flatter than the tracer density. If this region is weighted 
more in the virial integrals, then the average K zz required to 
produce the flattening is greater and so the virial sequences 
of an NFW halo fall below those of a logarithmic halo and 
also the observational data (cf. Figs. [5] and [6]). Thus, a log- 
arithmic potential with a finite core seems to be favoured 
over an NFW profile, at least as judged from the present 
constraints on the thick disc kinematics. 
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Figure 7. Virial sequences for the stellar halo in a spherical log- 
arithmic dark halo plus a short thin disc. The relative weighting 
of the disc and halo is controlled by fi = 0.7. Each virial sequence 
corresponds to fixing b = Rb/Rdm and varying d = -Rthin/^b 
from 0.09 to 0.12 (se e text). The models all lie above the pale blue 
ISmith et al.l l|2009al ) databox. Lower values of fj, yi eld somewhat 
higher pressure ratios, but still below the pale green iKeplev et al.l 
1120071 ) databox. The purple shaded zone shows models with con- 
stant anisotropy and spherical alignment, with the extreme ra- 
dial orbit model providing the diagonal edge. [Blue curves are 
for b = 7, purple for b = 8.6, green for b = 10.2 and orange for 
b = 12.] 

7 THE STELLAR HALO REVISITED 

In Section[5] we presented a simple model of the stellar halo, 
and showed that if the local measurements of the velocity 
dispersion using the SDSS subdwarfs (Smith et al. 2009a) 
are a good guide to the global kinematics, then this implies 
a moderately flattened potential with an axis ratio of g « 
0.85. If this is ascribed solely to the dark matter halo, then 
it corresponds to a flattening in the dark matter density 
contours of « 0.7. 

In Section [6l we showed that the global kinematics of 
the thick disc are consistent with motion in the combined 
potential of a spherical dark halo together with a thin disc. 
Here, we will ask: is it possible that the dark halo is spheri- 
cal, and that the flattening in the potential measured from 
the kinematics of halo stars is due entirely to the presence of 
the thin disc? If the dark halo is round, then the flattening 
of the combined disc and halo equipotentials changes from 
« 0.5 close to the Galactic centre to « 1 at Galactocentric 
distances beyond 10 kpc. 

The gravity field is now provided by the thin disc po- 
tential with short scalelength ~ 2.6 kpc, together with a 
completely spherical logarithmic dark matter potential as 
in eq l)47|l . The stellar halo density is given by the broken 
power-law of eq (|39[1 with break radius Rb- The two length- 
scales involved are 

b =RZ d= ~RV (58) 

We compute virial sequences for b in the range 7 to 14 and d 
in the range 0.09 to 0.12, which are suggested by the roughly 
known values of the scalelengths (that is, i? t hm = 2.6 ± 0.4 
kpc, R h w 2 - 4 kpc and Rb « 27 kpc) 

Fig. shows the pressure ratio Hrr/H zz as a function 
of q, the flattening of the stellar halo. The light blue and 
green boxes indicate the observational regimes, depending 
on whether the pressure ratio is estimated from the data 
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of ISmith et~ai1 ( |2009ah or of iKeplev et alj (|2007f ) assuming 
cylindrical alignment. In each c ase, the const r aints o n flat- 
tening come from the work of iDeason et all J2011bh . The 
virial sequenc es lie below t he Keplev et al.l (|2007h databox 
and above the lSmith et alj |2009al ) databox. Accordingly, it 
seems reasonable to argue that the underlying assumption 
of a spherical dark halo may hold good, as the models are 
not too far off the uncertain data. The purple zone covers 
the parts of the (Hrr/Uzz, q) plane for which the velocity 
ellipsoid is spherically aligned and the anisotropy is constant 
using eq (|27|l . The upper diagonal line corresponds to the 
purely radial orbit model. Hence, the models that do repro- 
duce the observations are radially biased with an anisotropy 
parameter f3 ~ 0.9. This makes good sense, as the under- 
lying dark halo is nearly spherical at large radii where the 
disc becomes unimportant, in which case the flattening can 
only be produced by substantial radial anisotropy. However, 
such radial anisotropy is much more extreme than implied 
by Kepley et al.'s (2007) or Smith et al's (2009a) data, which 
have P w 0.6. 



8 CONCLUSIONS 

The tensor virial theorem can provide a powerful technique 
for extracting gross structural properties, such as flatten- 
ings, from datasets. It has been widely used in galactic 
astronomy, particularly for self-gravitating systems - such 
as stars in elliptical galaxies — moving in the potential 
generated by their density field. The literature on appli- 
cations to tracer populations, in which the stars are pri- 
marily moving in the gravity field provided by dark matter 
or other stellar components, is quite spars e dWhitd fl989t 
ISommer-Larsen fe Christensen|[l989l ; iHelmillioOSl '). 

We have provided the first systematic exploitation of 
the tensor virial theorem for tracer stellar populations, and 
shown how it can be used to yield constraints on global kine- 
matics and flattening. This technique offers considerable ad- 
vantages over, for example, the Jeans equations. The latter 
depend on the derivatives of densities and velocity disper- 
sions, and are necessarily subject to uncertainties with noisy 
data. The tensor virial theorem deals with integrated bulk 
quantities, which are less susceptible to observational errors 
or systematic problems such as substructure in the data. 

An obvious application is to Population II stars mov- 
ing in the dark halo of our Galaxy. For spheroidal tracer 
populations with flattening q moving in spherical power-law 
potentials, there are analytic results that prescribe the ratio 
of globally averaged equatorial to vertical velocity disper- 
sion. Thanks to the Flattening Theorem, the results can be 
straightforwardly adapted to spheroidal dark halos with axis 
ratio g in the equipotentials. The virial quantities depend 
only on the ratio q/g in this instance. 

The main astrophysical results of the paper are three- 
fold. First, we have shown how to interpret the stellar kine- 
matical measurements of halo stars in terms of the global 
flattening of the dark halo. In particular, if the Smith et 
al. (2009a) values for the velocity dispersion of the halo 
subdwarfs are accepted, then the dataset is consistent with 
oblate dark halo models with a flattening in the potential 
of g ~ 0.85, corresponding to a flattening in the dark mat- 
ter density contours of « 0.7; this is driven by the radial 



anisotropy that Smith et al. found. If however, the Kepley 
et al. (2007) values are preferred, then the dark halo may 
well be close to spherical. We have to make the assumption 
that the local kinematics are proxies for the global kine- 
matics, which is assuredly a weakness in such analyses at 
present. However, the GAIA satellite, which is scheduled to 
fly in 2013, will for the first time provide us with detailed in- 
formation on the three-dimensional kinematics of halo stars, 
and so there is hope for real advances using this technique 
in the near future. 

Secondly, we have demonstrated that the dual halo 
structure as proposed by Carollo et al. (2010) appears to 
be incorrect. It is not possible to support an inner halo and 
an outer halo with the claimed velocity dispersions in the 
same gravitational potential: the model is inconsistent with 
the requirements of the virial theorem. A powerful aspect 
of virial methods is the fact that multiple populations re- 
siding in the same potential can be used to cross-constrain 
the dark hark halo mass and shape, without the need of 
building detailed models to reproduce the photometry and 
kinematics. There is considerable scope for extending our 
virial techniques to study multiple populations moving in 
the same gravity field, a topic th at is already important in 
studies of the dwarf sphero idals (|Amorisco fe Evans Il201ll ; 
IWalker fe Penarrubia fl201ll ). 

Thirdly, we have carried out an investigation into the 
kinematics of the thick disc. We find that the reported thick 
disc kinematics are consistent with a thin disc with scale- 
length and an isothermal dark halo. Specifically, we find that 
the thick disc kinematics and shape is consistent with a 
potential generated by a thin disc with scalelength « 2.6 
kpc, and surface-density normalisation of 56 ± 6M0pc~ 2 
(|Holmberg fc Flvnnl [2004 ). Although NFW potentials can 
be used to fit the rotation curve, they seemingly fail to repro- 
duce the virial ratios expected for the thick disc, as judged 
from the observational data. 

The near future sees a number of ambitious radial veloc- 
ity surveys, such as ESO-GAIA and Hermes. These will com- 
plement the abundant positional and proper motion data on 
stars in the Galaxy that will be delivered by the GAIA satel- 
lite. Given the quality and quantity of the datasets, it is an 
attractive proposition to look at virial quantities to gain in- 
sight into the Galaxy's potential, perhaps as a prelude to 
detailed phase space modelling. We anticipate that the ten- 
sor virial theorem will become an important tool to make 
sense of the very large datasets on stellar positions and ve- 
locities that will become available in the future. 
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